On effects of regular S=l dilution of S=l/2 antiferromagnetic 
Heisenberg chains by a quantum Monte Carlo simulation 



in 
o 
o 

(N 

Of 
< 



Fengping Jin, Zhaoxin Xu*, Heping Ying and Bo Zheng 
Zhejiang Institute of Modern Physics, 
Zhejiang University, 
Hangzhou 310027, P.R.China 

(Dated: February 2, 2008) 

The effects of regular S—l dilution of S=l/2 isotropic antiferromagnetic chain are investigated 
by the quantum Monte Carlo loop/cluster algorithm. Our numerical results show that there are 
two kinds of ground-state phases which alternate with the variation of S 1 = 1 concentration. When 
the effective spin of a unit cell is half-integer, the ground state is ferrimagnetic with gapless energy 
spectrum and the magnetism becomes weaker with decreasing of the S 1 concentration p = 1/M. 
While it is integer, a non-magnetic ground state with gaped spectrum emerges and the gap gradually 
becomes narrowed as fitted by a relation of A « 1.25-yp- 
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PACS numbers: PACS: 75.10.Jm, 75.40Cx 



I. INTRODUCTION 



The effects of substitutions of magnetic im- 
purities on the antiferromagnetic spin chain 
have attracted great interests in the past 
decade. It has been shown theoretically that 
the ground state properties vary with different 
dilution cases. For the random substitutions, 
the most interested case is of the S=l/2 impu- 
rities in Haldane chain pj. For example, inelas- 
tic neutron scattering experiment on the com- 
pound Y 2 BaNi0 5 substituting Ca 2+ for Ni 2 + 
[2] show a substantial increase of the spectral 
function below the Haldane gap to indicate the 
creation of states below the energy of the spin 
gap. This effects are also studied by numerical 
works by S. Wessel 3]. For regular substitu- 
tions, these systems are the mixed-spin chains 
which have been extensively studied by many 
authors in the past a few years. Analytical 
methods of non-linear sigma model, mean field 
theory and spin- wave method [IjOjlIlJ as well 
as numerical works by density matrix renor- 
malization group @ and quantum Monte Carlo 
[U have been applied extensively for such 
systems. So far, it is well known that the topol- 
ogy of spin arrangements in the mixed chains 
plays an essential role on the ground state prop- 
erties and thermodynamics in the mixed-spin 
systems. 

Experimently, many Quasi-ID mixed-spin 
materials have been synthesized in the past 
two decades, such as ACu(pba)(H 2 0)3 ■ n(H 2 



O) and ACu(pbaOH)(H 2 0) 3 • nH 2 (where 
pba= 1 ,3-propylencbis(oxamato) ,pbaOH=2- 
hydroxo-l,3-propylenebis and A= 

Ni,Fe,Co,Mn,Zn). These materials contains 
two different transition metal ions per unit 
cell, and their pro per ties were studied as fer- 
rimagnetic chains [U El El El El EH The 
experiment results imply that the magnetic 
properties of the mixed-spin compounds can 
all be described by a Heisenberg model with 
nearest-neighbor antiferromagnetic coupling 
as 
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H — JiSi ■ Si + \ 



(1) 



where Si denotes a spin-5 moment at site i, 
N is system size and Ji > 0. T. Fukui and 
N. Kawakami Q have studied spin chain com- 
posed by a periodic array of impurities S 1 em- 
bedded in the host S 2 ^ S 1 spin chain with the 
period M, i.e., 



S 1 ® S 2 (8) S 2 
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(2) 

The dilutions of the model denoted by the 
impurity concentration p has two limits: (i) 
p = 0, the undoped pure antiferromagnetic S 2 
chain, it has a non-magnetic ground state; (u) 
p = 0.5, the alternating spin chain of 5* 1 and 
S 2 . According to Marshall theorem and Lieb- 
Schultz-Mattis (LSM) theorem [T3|, ground 
state of the doped cases are specified by the 
spin quantum number S = 0(\S± — 5* 2 |iV/A/) 
for M = odd or even, it is either a spin sin- 
glet or ferrimagnetic. If the effective spin in a 
unit composed of M spins S e ff is half-integer, 
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so the system has a gapless energy spectrum. 
But when S e ff is integer, LSM theorem fails 
to predict the energy spectrum to be gaped or 
gapless. By applying non-linear a model, it is 
found that the system has an energy gap when 
the S e ff is integer 0,0]. But details of ground 
state properties and thermodynamics can not 
be given by non-linear a model analyses. 

The authors of present paper have recently 
studied the model (1) with the case of S 1 = 1/2 
and S = 1 by applying quantum Monte Carlo 
simulations [18j, where the numerical results 
reveal different non-trivial magnetic properties 
happened between two kinds of diluting cases, 
i.e. for odd S 2 — 1 spins in a unit, system 
has magnetic ground state and it shows ferri- 
magnetic features; while for even S 2 — 1 spins 
in a unit, systems behave non-magnetic ground 
states with antiferromagnetic-like features. For 
both the odd — even cases, the ground states 
are gapless steadily. And the system gradually 
transits from the ferrimagnetic ground state of 
the alternating S 1 -S 2 chain to the disordered 
ground state of pure S = 1 chain in two dif- 
ferent tendencies. In this Letter, we study 
an opposite case with 5* 1 = 1 and S 2 = 1/2. 
Previous analytical work predicted that if odd 
S 2 = 1/2 spins in a unit, the effective spin S e /f 
is half-integer, the ground state is ferrimagnetic 
with a gapless energy spectrum; while if even 
S 2 — 1/2 in a unit, S e ft is integer, the ground 
state is non-magnetic and the system has an 
energy gap. Our numerical study will focus on 
how the ground state properties depend on the 
concentration and the finite temperature mag- 
netic properties evolute with decreasing of the 
S 1 — 1 concentration p. 



II. CALCULATION AND RESULTS 

We use the efficient continuous imaginary 
time version of loop cluster algorithm to per- 
form the quantum Monte Carlo simulation |l9j , 
which has been successfully ap plied for the 
other mixed-spin chains [T^, |2(j . We confine 
our calculation to isotropic antiferromagnetic 
coupling cases, i.e. J; = J > in equation 
iJTJl. and the positions of spin Si = h and 
($2 = 1 are arranged as represented in equation 
(J2J) with M taking the values from 2 to 1 1 . We 
carry out 10 5 Monte Carlo steps for measuring 
physical quantities after 10 3 Monte Carlo steps 
for the thermalization. In order to clearly ex- 
plore the ground state properties, the simula- 
tions are performed at the very low tempera- 



ture (3 = 1/T = 200 for system sizes L > 200 
in condition of even number of unit. The phys- 
ical quantities we measure are the ground state 
energy Eg, the uniform magnetic susceptibil- 
ity Xu and staggered susceptibility \s by using 
the improved estimators in the loop cluster al- 
gorithm, e.g. , 



< X > 



= JL( 
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(3) 



(4) 
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where Wt (c) is winding number of cluster c, and 
|C| is the cluster size. The magnetization and 
staggered magnetization are estimated by 



< M > -- 



3(£^ 



) 2 ) 

/ MC 



and 



<M 2 >= (3(£(-irSf) 



2 > 



MC 



(5) 



(G) 



The energy gap A is also estimated in the way 
given by Todo, [2l| 



A = lim — -. 

L~,co i T ,Q\L) 



(7) 



where £ T) o is the correlation length in the imag- 
inary time direction. 

The results for magnetizations and uniform 
susceptibility are plotted in Fig. Q]and Fig. [5] 
We find that the magnetic properties are ap- 
parently different for two cases of M = odd and 
M = even. When M = even, the magnetiza- 
tion is finite and approaches zero linearly with 
decreasing of p. While M = odd, the magne- 
tization remains almost at zero value. On the 
other hand, it can be observed from our results 
that the uniform susceptibilities Xu is finite for 
M = even, but it vanishes when M = odd. 
Thus there is magnetic long-range order (LRO) 
in the ground state when M = even, but the 
order is absent when M = even. 

We further estimate the staggered magneti- 
zation and its susceptibility as a function of 
concentration shown in Fig. |3| and Fig. 0] re- 
spectively. The two observables are both finite 
for M = odd and M = even cases, but the data 
for the cases of M = even have much stronger 
values than the cases of M = odd. 

In order to confirm the results observed 
above, we begin to investigate the finite tem- 
perature uniform magnetic susceptibility. As 



3 



1 0.3 
bo 



FIG. 1: The magnetization versus p — 1/M, i.e. 
the diluting concentration of S 1 . The filled squares 
present the cases of M = even and the empty 
squares for M = odd. 



FIG. 3: The staggered magnetization versus p = 
1/M. The filled triangles are for the cases of M = 
even and empty ones for M = odd. 



FIG. 2: The magnetic susceptibility versus p — 
1/M. The filled diamonds are for the cases of M = 
even and the empty diamonds for M — odd. 



FIG. 4: The staggered susceptibility versus p = 
1/M. The filled triangles are for the cases of M — 
even and empty ones for M = odd. 



displayed in Fig. [5J one can easily find that \u 
diverges when the temperature T — 1/(3 goes 
to zero in the cases M = even. This is the typ- 
ical behavior of a system with magnetic LRO. 
In the cases M = odd, all the \u approach zero 
when T — > 0, a remarkable evidence to reveal 
the existence of the energy gap. 

Up to now, our results verify numerically 
that there are magnetic LRO and antiferro- 
magnetic LRO in the ground states when M = 
even. They clearly show that the ground states 
are ferrimagnetic in such cases. While for M = 
odd, there should exist of spin liquid phases 
denoted by the vanish of the magnetizations. 
Consequently we believe our numerical results 
consist correctly with the previous analytical 
predictions. More important, one can easily 
see that the magnetism decreases with decreas- 
ing of impurity concentration in the case of 
M = even. But there is not notable change 
of the magnetic properties when the S* 1 = 1 
concentration decreases as M — odd. 

Next, we consider the feature of the en- 



ergy gap A on different regular dilutions. Not 
surperised for us, the energy gap is closed when 
M is even and it opens again while M is odd 
as shown in Fig. |SJ These results is consist 
with the prediction by non-linear a model and 
LSM theorem 0, ■ It is interesting that the 
energy gap A tends to be narrow as decreasing 
of S 1 = 1 concentration when M — odd. We 
confirm such behavior by fitting A to the curve 
of 1.25-y/p as one can see in Fig. El 

Moreover, we show the finite-size effect of A 
results for several cases with L increasing in 
Fig. In our estimations, although the gaps 
are not exact closed for M = even due to the 
finite-size simulations, we find the data of the 
gaps decrease fast than L~ l , so it is obvious 
that the gaps will trend to zero as L — > oo. For 
the cases M = odd, where the gap opens all 
the time, there is almost no finite-size effect. 

In order to identify the ground state phases, 
we calculate the valence-bond-solid (VBS) 
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FIG. 5: The uniform magnetic susceptibility \ u 
versus temperature. The empty symbols are for 
the cases of M = even and filled ones for M = odd. 
The star present the result of undoped S = 1/2 
antiferromagnetic Heisenberg chain. 



FIG. 7: The finite-size effect of energy gaps. The 
two upper lines of filled symbols are for cases of 
M = odd whose values are presented by the y- 
axis at left-hand side, and the two lower lines of 
empty symbols are for M = even whose values are 
presented right-hand side y-axis. 
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FIG. 6: The energy gap versus a, i.e. the diluting 
concentration of S . The filled circles present the 
cases of M = odd and empty ones present the cases 
of M = even, the dashed lines are guide for eyes. 
The thin line is the fitting curve A = 1.25^^—0.1, 
the correlation coefficient of this curve to the empty 
circles is 0.999. 



order parameter 
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exp[*-5>J]>, 



(8) 



According to the LSM theorem, z vanishes in 
the gapless phase as system size N — ► oo. On 
the other hand, one expects that z varies in 
between ±1 but z ^ in a given gaped phase. 
In exact VBS states, z = ±1 |2^|. Our calcu- 
lations are plotted in Fig. |H1 It is clear that 
z « — 1 for all cases of M = odd to present the 
system located in a VBS phase; while z f» 0, 
it reveals the gapless energy spectrum for all 
M = even cases. 

Especially, all these ground state phases can 
be understood under the scenario of VBS pic- 
ture. In VBS picture, each impurity S 1 = 1 can 



FIG. 8: The VBS order parameter z versus 1/M, 
i.e. the concentration of S 1 = 1. The filled trian- 
gles are for the cases of M = odd and empty ones 
are for M — even. 



be regarded as two spin-1/2 in a triplet state, 
these two spin-1/2 can form singlet with their 
nearest neighbor S=l/2 spin due to the anti- 
ferromagnetic coupling. When M = odd, each 
unit have even number of S=l/2 host spins, 
so they can fall into singlets with their nearest 
neighbors including the two spin-1/2 of S = 1 
to induce the VBS order as seen in Fig. |5| (a). 
As a result, the system now shows a gaped en- 
ergy spectrum. But for M = even, odd number 
of spin-1/2 exist in a unit and there will be an 
active spin which is not used to form singlet as 
shown in Fig. HJb), thus there is no VBS or- 
der and the system emerges no spin gap. Our 
results of the VBS order parameter z clearly 
verify this picture, z ~ — 1 when M = odd and 
z w when M = even as shown in Fig. [S] At 
last, we note that the VBS phase is stable with 
the variation of S 1 — 1 concentration p when 
M = odd. 
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FIG. 9: The illustrations of VBS picture. The 
dashed ellipses present the spin 5 = 1, filled circles 
present the S — 1/2 and empty circles present the 
active 5=1/2 spins, (a) for M = odd and (b) for 
M — even. The empty circles in (b) present the 
active spin-1/2. 
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zero 


large 


small 


large 


<M S > 


finite 


finite 


finite 


finite 




small 


large 


large 


large 


A 


gaped 


gapless 


gapless 


gapless 


z 


-1.0 


0.0 


0.0 


0.0 



TABLE I: Comparison of ground state properties 
of two model, where S e ff is effective spin in a unit, 
<M> is magnetization, <M S > is staggered mag- 
netization, Xu is uniform susceptibility, Xs is stag- 
gered susceptibility, A is energy gap and z is VBS 
order parameter. 



III. DISCUSSION AND CONCLUSION 

Our Monte Carlo study verifies that two 
brunches of different magnetic behaviors 
emerge in cases of regular S=l diluted S=l/2 
host chains. According Marshall theorem, the 
cases with M — even have the ferrimagnetic 
ground states which can be specified by quan- 
tum number Stotai — \Si — S2\N/M, so the 
magnetization per site is finite and it decreases 
linearly as a function of p to M. = 0, the case of 
the pure S=l/2 antiferromagnetic Heisenberg 
chain. This feature can be easily observed from 
our results in Fig. ^ an d Fig- 13 When M is 
odd, the ground state is singlet with Stotai = 0, 
thus the magnetization per site keeps zero and 
this is a non-magnetic state. As observed in 
our simulations there is no notable variations 
of the ground state magnetic properties in the 
cases of M = odd. 

To compare our ground state results of the 
model in this Letter for S 1 — 1 and S 2 = 1/2 
(system T) with the one we studied previously 
when S" 1 = 1/2 and S 2 = 1 (system II), we 
collect the main points of the numerical cal- 
culations in Table I. One can easily see that 
both systems behave with two kinds of dif- 
ferent ground state phases, magnetic or non- 
magnetic, respectively. If M — even, the 
ground states are ferrimagnetic for both system 
I and system II, and their magnetizations and 
staggered magnetizations are all finite and de- 
crease linearly with decreasing of impurity con- 
centration. However, for the cases of M = odd, 
there appears VBS order in system I which is 
gaped, but the order is absent in system II 
where the spin arrangements can not induce 
such order, so the gap is constantly closed. 
This feature reveals that this topological or- 



der plays an important role to the behavior of 
the energy gap in the mixed-spin system. We 
believe that the fitted relation of Aw 1.25^/p, 
to denote the energy gap as function of S 1 = 1 
concentration, provides a good stuff to study 
how the topological order affects the energy gap 
in the mixed-spin systems. 

In conclusion, we have studied the ground 
state and finite temperature magnetic proper- 
ties of the regular S 1 = 1 diluted in S 2 — 
1/2 antiferromagnetic chain. Our calculations 
show that there exist different phases in the 
ground state as a function of S 1 concentration. 
When there is one S 1 impurity and odd num- 
ber of host S 2 spins in a unit cell, the ground 
states are ferrimagnetic and the system has a 
gapless energy spectrum. The ferrimagnetism 
becomes weaker as the impurity concentration 
reduced. While for one S 1 and even number 
of S* 2 in one unit cell, the ground state is a 
VBS phase where there is a gaped energy spec- 
trum and the energy gap gradually approaches 
to zero with decreasing the concentration p. An 
interesting observation is that the behavior of 
the energy gap can be numerically well fitted 
by A w 1.25y/p. Further analytical work, for 
exmaple using the mean-field theory 0, is re- 
quired to explain why such dependence of the 
energy gap exist in VBS phases. 
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